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Abstract 

A method is developed to construct the solutions of one and many variable, 

linear differential equations of arbitrary order. Using this, the A^-particle 

Sutherland model, with pair-wise inverse sine-square interactions among the 

particles, is shown to be equivalent to free particles on a circle. Applicability 

of our method to many other few and many-body problems is also illustrated. 
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The Calogero-Sutherland model (CSM) |l|,^ and its generalizations are extremely 
well-studied quantum mechanical systems, having relevance to diverse branches of physics 
0, such as the universal conductance fluctuations in mesoscopic systems 0, quantum Hall 
effect 10], wave propagation in stratifled flelds P, random matrix theory P,|5|J9|JT0[] , fractional 
statistics |TT]], gravity and gauge theories [0. In particular, the Sutherland model [Q, 



describing A^-identical particles having pair-wise inverse distance square interaction on a 
circle has received wider attention in the context of questions related to statistics The 



spectrum of this exactly solvable model can be interpreted as arising from a set of free 
quasi-particles satisfying the generalized exclusion principle [|15| . 



Some time ago, the present authors have shown that, the CSM is equivalent to a set 
of free harmonic oscillators [|1^]. This mapping throws considerable light into the algebraic 



structure of this correlated system ||T^,|T8|. Later, the same model, without the harmonic 
confinement, was made unitarily equivalent to free particle [jl9[. However, the connection 



between the Sutherland model and free particles have so far remained an open problem. In 
this paper, we establish the equivalence of the Sutherland model to free particles on a circle, 
for arbitrary number of particles. 

This paper is organized as follows, (i) A general method is developed to solve any 
linear differential equation of arbitrary order, (ii) a few examples are discussed to show how 
the present technique yields the known, as well as, new results, (iii) the usefulness of this 
technique for solving the Schrodinger equation with complicated potentials is pointed out, 
(iv) the mapping between the A^^i Calogero model and free harmonic oscillators is rederived 
using the present method, and, finally, (v) we show, the equivalence of the Sutherland model 
with free particles on a circle. 

Consider the following differential equation, 

{F{D) + P{x,d/dx))y{x) = , (1) 

where, D = x^, F{D) = J2n='^oD(^nD"', is a diagonal operator in the space of monomials 
spanned by and a„'s are some parameters. P{x,d/dx) can be an arbitrary polynomial 



function of x and We prove the following ansatz, 



lm=0 



FiD) 



P{x, d/dx) 



X 



(2) 



is a solution of the above equation, provided, F{D)x'^ = and the coefficient of x^ in 
y{x) — Cxx^ is zero (no summation over A); here, C\ is a constant. The case, when the 
equation F{D)x^ = does not have distinct roots is not considered here and will be treated 
separately elsewhere. 

Substituting Eq. (^, modulo Cx, in Eq. (|I]), 
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{F{D) + P{x,d/dx)){Y.{-l] 
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(3) 



Eq. (^, which connects the solution of a given differential equation to the monomials, 
can also be generalized to many- variables. In order to show that, this rather straightforward 
procedure indeed yields non-trivial results, we explicitly work out a few examples below and 
then proceed to prove the equivalence of the Sutherland model to free particles on a circle. 

Consider the Hermite differential equation, which arises in the context of quantum har- 
monic oscillator, 



(4) 



Here, F{D) = D — n and F{D)x'^ = yields X = n. Hence, 



oo 



-l/2){(f/dx'^ 



Using, [D , {(f/dx"^)] = —2{d'^/dx'^), it is easy to see that 
1 



X 



-l/2){dydx^' 



(D-n) 



X 



m 1 

-i/2r{d'/dx'rx{j-^^x^ 



and 



(5) 



Hn{x) = CnY. i-i/^r—idVdx^rx 

m=0 



_1 d 

C„e 4^2;" 



(6) 



this is a well-known result. Similar expression also holds for the Lagurre polynomials which 



matches with the one found in |2^. In order to make an important remark, we list the 
solutions of some frequently encountered differential equations in various branches of physics 
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Legendre polynomials 



Associated Legendre polynomials 



Bessel functions 



Generalized Bessel functions 



Gegenbauer polynomials 



4 



Hypergeometric functions 



y± 



where, X± is either a or (3 and A = x-^ + 'j-^. All the above series solutions have descending 



powers of x. In order to get the series in the ascending powers, one has to replace x by Mn the 
original differential equation and generate the solutions via Eq. (^. However, the number 
of solutions will remain the same. One can also generate the series solutions by multiplying 
the original differential equations with x^, and then, rewriting x^-A- = D(D — 1) = F{D). 



Meijer's G- Function 

{oo 
E 
A:=0 

Neumann's polynomials 



(-1) 



p—m—l 



xY[{D — Gj 



X 



OJx) 



E(-i) 



.r=0 



■X 



j| (a;cos^(n7r/2) 
+?T, sin^(r;,7r/2)) ; 



the cases of Struve, Lomel, Anger and Weber functions are identical to the above one. Jacobi, 
Schlafli, Whittaker, Chebyshev and some other polynomials were not given here since the 
list is rather lengthy. Now, the remark follows: Gx becomes independent of the roots of the 
equation F{D)x^ = 0, only when F{D) is linear in D. 

The solution for the following equation with periodic potential, 



dx 



Uj y 

-j-^ + a cos{x)y = 



(7) 



can be found, after multiplying Eq. by and rewriting as {D — 1)D, to be 



y{x) 



°° ^-a)'" (-1) 



E 



m,{ni}=0 



n 



n 



{2[m + \/2-r + Y.T=i"'n,])\ 



X a;2(m+^-^n,+A/2) ^ (g) 



where, A = or 1. In the same manner, one can write down the solutions for the Mathieu's 



equation as well ^ 



Using this method, one can, in principle write down the ground and first excited states 
of the one variable Schrodinger equation. We first illustrate the procedure for harmonic 
oscillator, {h = m = u = 1), 



^ + (2E„ - X^Yj = 



(9) 



before proceeding to the other non-trivial examples. Multiplying by x^, the above can be 
written as 



[{D - 1)D + x'{2En - x^)) ^„ = 



(10) 



For n = and 1, D{D — = 0. Using Eq. (j^), the solution for n = is, 

1 
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{D - 1)D 



{x'{2Eo-x')) 



X 



[2Eo] 2 , (2! + [2Eofj , 



4! 



X — 



(4! + (2!)2[2Eo] + (2!)[2Eo]^ B 



2!6! 
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Note that, ipo is an expansion in the powers of x, whose coefficients are polynomials in Eq. 
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Only when Eq = 1/2, the series can be written in the closed form Coe~2^ , which satisfies all 
properties required for a quantum mechanical ground-state. Similarly, the first excited state 
can also be found by applying Gi on x. To find the n-th excited state, one can differentiate 
the Schrodinger equation n times, multiply it by and use x"'-^ = 1^=0 {D — I) = F{D). 
In the case of complicated potentials, this method can be potentially very useful if applied 
in conjunction with numerical algorithms. The exact ground-state wavefunction for the 
Schrodinger equation with anharmonic potential, V{x) = x^/2 + cx'^/2, is, 

1 



k m=0 



-{x\2Eo - x^ - cx^)) 



.1 



(12) 



iD-l)D 

One can numerically tune Eq to obtain an appropriate ipQ. 

In the following, we apply the above technique to the differential equations involving 
many variables. Consider, the A^^i Calogero-Sutherland model, which was made equivalent 



to free harmonic oscillators |16], 
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N 



d' 



N 
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1 



V 



En {i'0(j)n 







(13) 



This can be brought to the following form, after the removal of the ground-state wavefunc- 
tion, ^0 = exp{-iEia;f}ni^j[ki - Xjl": 

d 



Y^Xi-^ + Eo- En- A 







(14) 



where, Eq = + ^N{N — l)a is the ground-state energy and A = [\Y.i Tpr + 



1 a 



*7^i (Xi—Xj) dXi 



\. Without loss of generality, we have quantized the above system as bosons. 



Rewriting Eq. (|T^ as. 



d 



^ - + n + - -Sn - ^ 



(15) 



the solution can be written by using Eq. (H): 

1 



E(-i; 



m=0 



n + Eo-En-A) 



Sn{{Xi}) 



(16) 



here, iJ^iXi-^ — n)Sn{{xi}) = and S'„({xi})'s are any homogeneous function of degree n. 
In order to avoid the possibility of having singular solutions when the inverse of {Yl,i ^«^~''^) 
acts on Sn{{xi\), we choose n + Eq — En = this yields the familiar energy spectrum of the 
Calogero-Sutherland model. Further, one has to choose Sn{{xi}) to be completely symmetric 
under the exchange of Xj's, such that, the action of A yields polynomial solutions which are 
normalizable with respect to ipo as the weight function |ll6| , |22[| . Similar to the case of the 
Hermite polynomial, one can easily prove that. 



0n = Cne-^^Sn{{Xi}) 



(17) 



In the following, we show that, the Sutherland model is equivalent to free particles on a 
circle. The Schrodinger equation is. 



N ^ 

dx'j 



1 



EAM{X^})=0 



(18) 



i=l 
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Choosing, Zj = e^'^*^^/'^ and writing ip\{{zi}) = Hi 



-(Af-l)/3/2 



Ui<j{zi - Zj)'^Jx{{z^}), the 



above equation becomes, 



Zi ~t~ Zj 



i<j 



{D,-D^)+Eo-EAJx{{z^})=0 



(19) 



where, A = = (^)'^a, = (^)'^o and = |(f )'/52iV(iV2 - 1), is the ground- 



state energy. Here, Jx{{zi}) is known as the Jack polynomials p3H26|l. Y^iD^ is a diagonal 
operator in the space spanned by the monomial symmetric functions, with eigenvalues 

E^i Af Rewriting (|g) in the form. 



(20) 



one can immediately show that. 



n=0 



E ^^(A - /^,) + ^0 + E A' - A 



xmA({2;i}) 



= CxGxmx{{zi}) 



(21) 



It is easy to check that, Gx maps the Sutherland model to free particles on a circle, i.e., 

{^PoGx)-'Hs{^oGx) - 



271 

T 



)'(ea^-Ea.^ + a) 



N 02 



(22) 



where, Hs is the Sutherland Hamiltonian and ipo is its ground-state wavefunction. For the 
sake of convenience, we define 

1 



S 



and Z^/35:^^(A-/^,) + ^o + E\'-A 



(23) 



The action of S on mx{{zi}) yields singularities, unless one chooses the coefficient of nix in 
Z mx{{zi\) to be zero; this condition yields the well-known eigenspectrum of the Sutherland 
model: 

Ex = E^ + E(\' + m+^- 2^]Ai) 
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Using the above, one can write down the Jack polynomials as 

1 



n=0 



xmx{{xi}) 



(24) 



For the sake of illustration, we give the following computation of J2 for two-particle case. 
Now, 1712 = + Z2 and it is easy to check that. 



Zm2 = l3AziZ2 = A(3mi2 

Sm2 = ^^(^^(4/5-1^) 
^"m2 = -2(/3)"mi2 for n > 1 . 



Substituting the above result in Eq. (piD, apart from C2, 

J2 = m2+ ff;(-ir(-2)(/3r) m,. 



\n=l 



m2 + 2f3[Y^{-f3r]m,. 



\n=0 



m2 + ^m,. 



(25) 



this is the desired result. Earlier, Lapointe et al. have obtained a Rodrigues-type formulae 
for the Jack polynomials Eq. (^) is a new formula for the Jack polynomials. In the 
following, we show that, the above method can be generalized to more general many-body 
systems. 
Consider, 



E «n(EA")+^M^A({^.})=fiA({^.}) 



(26) 



where, a„'s are some parameters, Di is same as the above, A is a function of Zi and d/dzi 
and B\{{zi}) is a source term. 

Case (i) When Bx{{zi}) = and Anix = exmx + J2fi<x C^xttIii] where, w-a's are the monomial 
symmetric functions [^] and eA and Cxfi are some constants, we get 

1 



r=0 



mxi{z,}) (27) 



with, E^=-oo«n(EiAr) + eA = 0. 

Case (u) When Bx{{zi}) = J^t^^xK^m^ and Am^ = E^^^a C'^a^^m^ where, and Cxf, are 
some constants, then, 

1 



r=0 



1 



5a({^.}) . (28) 



In conclusion, we have developed a general method to solve any linear differential equa- 
tion of arbitrary order. The major advantage of this technique lies in the fact that, in order 
to obtain the solutions, one does not have to solve any recursion relations or integral equa- 
tions. By applying this method, we obtained new formulae for many differential equations 
which are of interest to theoretical physicists. We also discussed how to obtain the ground 
and excited states for the Schrodinger equations with complicated potentials. In particu- 
lar, we have written down the exact ground-state wavefunction for anharmonic potential, 
V{x) = + as a function of energy parameter, which has to be tuned to obtain 

the normalizable solutions. This method treated both A^v-i Calogero-Sutherland and the 
Sutherland model with inverse sine-square interactions on the equal footing. The former one 
is equivalent to free harmonic oscillators and the later, to free particles on a circle. We also 
proposed a general scheme for generating more general symmetric functions. This scheme 
may find interesting applications in the context of obtaining the g-deformed versions of the 
known polynomials pT|j25|] . 

The authors acknowledge extremely useful discussions with Profs. V. Srinivasan and S. 
Chaturvedi. N.G thanks U.G.C (India) for the financial support. 
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